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0,1 We study the evolution of non-Gaussianity in multiple-field inflationary models, focusing 

I 1 on three fundamental questions: (a) How is the sign and peak magnitude of the non- 

O ' linearity parameter /nl related to generic features in the inflationary potential? (b) How 

sensitive is /nl to the process by which an adiabatic limit is reached, where the curvature 
C^ , perturbation becomes conserved? (c) For a given model, what is the appropriate tool — 

C^ ■ analytic or numerical - to calculate /nl a-t the adiabatic limit? We summarise recent 

results obtained by the authors and further elucidate them by considering an inflection 

point model. 
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1. Introduction 



An important feature of canonical single field models is that the curvature pertur- 
bation produced at horizon crossing is conserved, ^'^ with statistics indistinguishable 
from Gaussian.'^ In multiple field models, on the other hand, isocurvature modes 
may also be produced, and can subsequently source the evolution of the curvature 
KS ' perturbation as the field space path curves."*'^ The curvature perturbation and its 

j^ ■ statistics (such as the power spectrum and non-Gaussianity®'^) can therefore con- 

tinuously evolve after horizon crossing, leading to far richer behaviour than that 
allowed in single field models. Here we focus on non-Gaussianity, and in particular 
the nonlincarity parameter, /nl, though many of our general conclusions extend to 
other obscrvablcs as well. 

In canonical multi-field models, the non-Gaussianity present at horizon cross- 
ing is negligible.® In this setting, however, it is possible for /nl to evolve to large 
values. ^'^°'^^'^^ To calculate the observationally relevant value of non-Gaussianity, 
in principle one has to follow the evolution until the time of last scattering, where 
the Cosmic Microwave Background (CMB) was imprinted. Given our present igno- 
rance about the detailed physics of the early universe, this would not be possible in 
practice. In many models, however, a regime is reached during the evolution long 
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before this time - the so-called adiabatic limit - where all the isocurvature modes 
decay, and the curvature perturbation becomes conserved. 

There are different ways in which such a limit could be reached, and this has 
consequences both for the possible values of the observable parameters, such as /nl, 
and for the techniques which can be reliably employed to calculate observables, i.e. 
whether analytic methods will suffice, or numerical methods are necessary. A useful 
classification of models according to when the adiabatic limit is attained is: 

• Models in which an adiabatic limit is reached 'naturally' by convergence 
into a 'focusing region' of the potential. This can be further sub-divided into 
cases where the convergence occurs during slow-roll, and cases in which 
convergence occurs only after the slow-roll approximation fails. 

• Models where an adiabatic limit is reached abruptly due to an additional 
degree of freedom, such as a waterfall field being destabilised. 

• Models for which no focusing region in the inflationary potential exists, and 
an adiabatic limit can be reached only by embedding the inflationary model 
into a larger scenario, perhaps one which includes perturbative reheating. 

Recently an extended study of these possibilities was undertaken by the authors 
in the context of of multi-field models of inflation, capable of producing large non- 
Gaussianities.^^ Here we give a summary of those results and further elucidate them 
by considering an inflection point model, which is of the flrst type, and which can 
produce a large positive or negative value of /nl at the adiabatic limit. 

The structure of the paper is as follows. In §2 we give a summary of the back- 
ground theory which is used to formulate analytic expressions for observables, dis- 
cussed in §3, together with the conditions needed for the non-Gaussianity to be 
large at the adiabatic limit when it is reached naturally. In §4 we briefly discuss 
features in the potential which generate large transitory non-Gaussianities during 
the super-horizon evolution, which may be relevant for the final observable value at 
the adiabatic limit if this limit is reached abruptly. Finally, we discuss the useful- 
ness of our results in §5 when applied to specific models, and demonstrate this by 
considering a new example. We conclude in §6. 



2. Background 

We consider inflation driven by multiple canonical scalar flelds (j)i with i = 
1,2,..., M, self-interacting through a potential 1^(01, 02, ■•■)• Deflning W^i — 
dW/d4>i the scalar equations of motion are 

0, + 3H^, + W^^ = 0, (1) 
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where H is the Hubble rate, given by the associated Friedmann equation SMpjiJ^ = 
W + ^j 0I/2. We define slow-roll parameters as 

such that for inflation we require e < 1. The 'slow-roll limit' is given by e ^ 1, during 
which the fields' kinetic energy may be neglected, the decaying modes discarded, 
and the field equations well approximated by 

ZH^, + W^i = 0, 3Af 2jij2 ^ ^^ ^3^ 

Primordial cosmological perturbations are commonly characterised in terms of 
the curvature perturbation on uniform density spatial hypersurfaccs, denoted by 
Q. An important feature of C, is that for adiabatic perturbations it is conserved on 
large scales,^ at the linear order and even beyond. In multi-field models, on the other 
hand, C can evolve due to the presence of isocurvature modes, and this may result 
in the production of large non-Gaussianities. Deviation of the three-point function 
from zero is commonly measured in terms of the dimensionlcss parameter /nl 

_ 5 klklkl g^(fci,fc2,fc3) 

^"''^ Qkl + kl + kl An^Vl ' ^' 

where Vq is the power spectrum and Bq the bispectrum, given respectively by 

27r2 

(Ck.Ck.) = {2i^f5\]^^ + k2)-^Vciki) , (5) 

fci 

(Cki Ck. Cks) = (27r)3j3(j^^ ^ k^ ^ k3)Sc(fci, fc2, fca) . (6) 

A common technique for calculating (^ and its statistics, including /nl, is the 
5N formalism, "'^^'"'^^'"'^^ based on the separate universe approach to perturbation 
theory.^'^^ In this approach, spatial gradients are neglected on scales greater than 
the horizon size, and each spatial point is assumed to evolve as a separate FRW 
universe. In phase space, this can be represented by a bundle of trajectories, each 
evolving along an independent path from perturbatively different initial conditions. 
The variables which parametrise this phase space are the scalar fields {(j)i,(j>i}, as 
well as any radiation or matter species that may be present. The idea of associating 
inflationary perturbations with trajectories has a long history.^'^^'^'^® 

In this picture, choosing a different spatial slicing of the universe corresponds to 
taking a different cross section of the bundle in phase space. 

Choosing a flat initial slicing at i = t*, and a later uniform density (constant H) 
slicing at t = tc, then C on the final slicing can be equated with the difference in the 
number of e-folds, as measured along different trajectories in the bundle, ( = 5N. 

In general, as the bundle evolves in the field space, so will 6N and its statistics, 
requiring the dynamics to be followed indefinitely. If, however, the trajectories con- 
verge to a line parametrised by a single variable, i.e. the adiabatic limit, ^ becomes 
conserved. Moreover, in this limit each value of the Hubble rate corresponds to a 
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single combination of the field or fluid content, implying that derivatives of final 
quantities on c, with respect to changes in the initial conditions, tend to zero as the 
adiabatic limit is reached. In particular, one finds 90^/90^ -> 0. 

Restricting our attention to purely scalar field dynamics, a common example 
of how an adiabatic limit is reached is for the trajectories to evolve into a focus- 
ing region of the potential, such as a potential valley, possibly terminating in a 
minimum. For convergence into a valley, the mass-squared matrix associated with 
perturbations orthogonal to the direction to which the trajectories are converging 
should have large and positive eigenvalues. Taking the smallest eigenvalue to be of 
magnitude ^ m±, one typically expects a decay of the field derivatives d(j)1 / d(j)*i. at 
least as fast as e"^'"-^^^/^^^-' (see our recent work-'^'^ for a detailed discussion). As 
was mentioned in the introduction, however, an adiabatic limit could be reached in 
other ways. It may turn out that no focusing region in the potential is available, in 
which case the decay of the fields into radiation may need be considered in order 
for the model to make unambiguous predictions. Alternatively, an adiabatic limit 
could occur due to a sudden transition, such as a waterfall field being destabilised. 

In any case, before we can consider the fate of observables at the adiabatic 
limit, it is first necessary to have calculable expressions for these quantities, which 
could be obtained using the 5N formalism discussed above. We recall that during 
slow-roll infiation field velocities are functions of field positions. Taking this to 
be a good approximation at horizon crossing, the subsequent number of e-folds 
undergone by any 'separate universe' is then a function purely of the initial field 
values, iV((/)*, ..., 0^), even if the evolution subsequently evolves away from slow 
roll. Taking t* as a time shortly after observable scales left the horizon, therefore, 
a Taylor expansion 

C^SN^Y. ^^^^^t + ^ E N^^jS^*S^*j +■■■ (7) 

i ij 

can be made, where TV is the number of e-folds from * to c, a subscript i represents 
a derivative with respect to cj)* , and S4>1 are the field fluctuations on the flat hyper- 
surface at horizon crossing. Such a Taylor series allows the properties of a bundle 
of trajectories to be parametrised by just a few numbers, namely the derivatives 
of TV about some typical member of the bundle. Moreover, Eq. (7) allows various 
statistics to be estimated. In particular /nl is given by^^ 

3. Analytic Schemes and /nl at a 'Natural' Adiabatic Limit 

To analytically evaluate the non-linearity parameter /nl from Eq. (8) we must cal- 
culate N,i and N^ij at time t'', or when they become constant at the adiabatic limit. 
Currently, analytic calculation is only possible when the slow-roll equations of mo- 
tion, Eq. (3), are a good approximation. This means that if (^ is still evolving at the 
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end of inflation, we cannot analytically follow observable quantities, and numerical 
simulations will become essential^. Moreover, calculations require a special 'sepa- 
rable' form for the potential^. Vernizzi & Wands^ and later Battefeld & Easther^^ 
studied sum-separable models, W = 'Ylii^ii.i'i)^ deriving expressions for these co- 
efficients and for /nl- Similar techniques were used by Choi et al?^ for models of 
product-separable form W = IiiVi{(t)i), and recently Wang^^ generalised the study 
of sum-type potentials to those of the form W = (Yl,i'^i^4'i)Y^^^ where A is an 
arbitrary constant*^ A summary of analytic expressions for Al-field models of these 
forms is given in our paper. ^'^ 

The analytic formulae follow by using Eq. (3) to write N as an integral over one 
of the fields 0fc as iV = — /, . W/{M'^^^^^ ) d(j)k ■ In general, taking the derivative of 
this expression yields three contributions, namely initial and final boundary terms 
and a path term. For potentials for which analytic progress is possible, however, the 
path terms are either absent, or the integrals can be manipulated to make them so. 
For potentials with product-separable forms one finds 



N'r' - 



Vk 



M^Vl 



Vk 



Hk 



Ki^ 



(9) 



where the free index k labels the Ai ways of writing N^i, all of which will lead to the 
same result once the c-dependent terms are evaluated. The summation convention 
is not used anywhere in this paper. For potentials of generalised sum separable form 
one finds a similar expression 



N_,, 



V 



M^Vl 



M 

E 



Vk 



Ml,Vl 



(10) 



In both cases A^ ,,, follows by differentiation. 

The difficult step in deriving an analytic expression for N^i is the calculation of 
the c-dependent derivatives in Eqs. (9)-(10). A case where analytic progress is much 
easier occurs if the adiabatic limit is reached during slow- roll infiation. If this limit is 
reached by convergence into a valley, ^^ we expect the c-dependent derivatives to tend 
to zero at least as fast as the lightest isocurvature mode decays. ^^ If sufficient time 
is available for them to become negligible, this greatly simplifies the expressions for 
N^i. Indeed it is often possible to set the entire c-dependent boundary term to zero, 
and the expressions become dependent only on the values the fields took at horizon 
crossing. Where this simplification has been used in the literature, it has been 
referred to as the Horizon Crossing Approximation (HCA).^^'^^ Caution is needed, 
however, since it is possible that the coefficients VkjVl in front of the derivatives 
may diverge as the adiabatic limit approaches. This possibility is discussed at length 



^We arc not considering models such as the curvaton, where approximate analytic formula can be 
derived for regimes after the end of infiation by modelling the curvaton field as a fluid. 
''The product-separable ^"^ and the generalised sum separable potentials'^^ were explicitly only 
considered with two fields, but the results are easily generalised to an arbitrary number. -'^'^ 
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elsewhere, ^'^ and is only possible if the field, (t>k^ is completely orthogonal to the final 
adiabatic direction (the valley bottom), and T4 tends to a constant as VJ, tends to 
zero. In this case the c-dependent term will tend to an unknown constant (unless 
the full calculation of the derivatives is performed), rather than zero. In the case 
of sum-separable potentials, however, such a possibility can be avoided since we 
are free to reparametrise the potential and associate the problematic constant with 
another field. In the product-separable case we must simply pick the fcth version 
of Eq. (9), associated with a field (pk which is not orthogonal to the final adiabatic 
direction. Once this procedure is followed, the second derivatives of N follow by 
differentiation'^. The technical details are presented in our paper. ^'^ Here we instead 
illustrate this point with the help of some examples. 

Consider first a sum-separable potential of the form W = Vo{e^'^ + e'*'-^), which 
has a valley bottom defined by the line (p = x- Since neither field is orthogonal 
to this final direction, both V^ and V' remain non-zero as the adiabatic limit is 
asymptotically approached, and both fields continue to evolve. As the c-dependent 
derivatives vanish, therefore, so does the entire c-dependent boundary term, leading 
to expressions for N_i which depend only on initial conditions. Alternatively, con- 
sider a potential of the form W = Wo + ^m'j^cj)'^ + ^m'^x^ with m^ » m^. In this case 
the adiabatic limiting trajectory is the line defined by x = 0. As convergence to this 
trajectory occurs, V' vanishes and, if we define V^ — Wq + ^rn^^x^ > "^6 would arrive 
at an incorrect expression for the N^i by naively setting the c-dependent terms to 
zero. We are free, however, to define instead V^ = im^x^ and V^ = Wq + \m?,<j)^, 
and with this definition we can correctly set the c-dependent term to zero. 

An example of a product-separable potential with an convergent valley region 
is W^ = Wq{\ + g4>^) exp(— Ax^). In this case an adiabatic limiting trajectory occurs 
when = 0, where the velocity of the field tends to zero. Following the discussion 
above, we select the way of writing iV.i associated with the field which is still evolving 
as the valley bottom is reached, in this case x, and with this choice setting the c- 
dependent term to zero leads to an accurate asymptotic expression for N^i. 

In this way any model with a separable potential can be analysed, and one im- 
mediate question is of interest: Can ff^h be large after natural focusing? Discarding 
the c-dependent terms one can employ Eq. (10) together with Eq. (8) to obtain for 
the product-separable potentials the expression 

^/nl = 26*-,7*,. (11) 

where 0^ is the field still evolving at the adiabatic limit. Thus we find that mod- 
els with product-separable potentials lead to slow-roll suppressed values of non- 
Gaussianities if the adiabatic limit is reached during slow-roll inflation. 

The situation for sum separable potentials is more complicated, but a relatively 
simple picture emerges if we assume that Nj is much larger for one field, 4> say 

■^Formally, one should first calculate the expression for N^ij, and then allow derivatives of ij>i to 
tend to zero. These two operations commute provided (Vj/V)' is finite at the adiabatic limit. 



/NL = -7^A/pV/00— • (12) 
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(or for just a few fields as detailed by Kim et al}^) In this case V^/VL at horizon 
crossing is much greater than the analogous terms for the other fields, and one finds 

Thus in single field inflation /nl would be suppressed by the slow-roll parameter 
ry. However, in multiple-field models Vj/ /Vfj, need not be small even when Ty^^ is, if 
one or more of the other fields contribute significantly to the energy density. The 
condition for a large non-Gaussianity in this case is therefore that the mass-squared 
of (j) is much greater than magnitude of its potential at horizon crossing (in addition 
to the condition on V^/V^). We note that the sign of /nl is opposite to that of Vf'. 

4. The Magnitude and Sign of Transient Non-Gaussianities 

Our primary interest is in the final constant value of /nl after an adiabatic limit has 
been reached. However, 'transitory' large values of /nl may also be relevant, since 
it is possible for the infiationary dynamics to be interrupted, rapidly establishing an 
adiabatic limit. Examples may occur in models containing a waterfall field In that 
case, an evolving but large /nl value can be preserved at the adiabatic limit. The 
general conditions for such a large evolving value of /nl were given by Byrnes et 
al.,^'^ using the full analytic expressions for separable potentials. ^'^° In our recent 
work"'^'^ a different point of view was adopted. We studied the corresponding condi- 
tions which are required for a large /nl to be produced by features that commonly 
occur in multi-field potentials. We developed an approximation, based on intuition 
from the phase-space picture of inflationary trajectories, which allows simple scal- 
ings to be derived for the peak magnitudes of the transitory /nl as well as their 
expected signs. Space does not permit us to detail the entire derivation, so here we 
briefly summarise the conclusions found. 

To date models studied in the literature, which are capable of producing large 
transitory /nl, possess two broad features in their potentials: a ridge or a valley, 
with a large /nl produced as the bundle falls from a ridge or begins its turn into 
the bottom of a valley. A common feature present in both these cases is the rota- 
tion of the bundle of trajectories which results in sourcing the evolution of C from 
isocurvaturc modes. Moreover, both cases lead to a highly non-linear dependence 
of N on initial conditions in the early stages of the turn. This is because one side of 
the bundle of trajectories begins to turn before the other, leading to a temporary 
asymmetry in the dynamics. In general this non-linear dependence is encoded in 
the SN formalism through large values of the second derivatives N^ij. 

To demonstrate this behaviour more concretely, we specialised to two-fleld (0, x) 
potentials and considered generic ridge and valley potential forms obtained by per- 
turbative expansions about some position along a separatrix or the bottom of a 
valley, respectively. 

For the ridge, the relevent potential takes the form^'^ W — Wq + g(f) — ^m'^x^i 
where Wo, g and m-^ are constants, and x — defines the ridge. We assumed that 
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the initial field position was sufficiently close to x = so that the evolution was 
initially almost entirely in the 4> direction. Employing the 5N formalism we then 
showed that the evolution leads to a negative peak in /nl with its maximum value 
scaling as^^ 

f I n-^^i/^-^pi (^■i\ 

A* 

This inverse scaling with x* ^^^ explicitly verified using a number of potentials 
containing ridge features. ^^ Two important points are worth noting here. Firstly 
the sign of /nl resulting from a ridge feature is negative and secondly, a large value 
of /nl requires a high degree of fine tuning of the form x* ~^ 0. 

In general the valley evolution is more complicated. However, a similar picture 
emerges when a perturbative expansion of the form W — Wq + ^to^0^ + ^m^x^j 
is considered where Wq, ?ti0 and m^ are constants. ^^ In this case we assumed 
m0 ^ m^, so that the initial evolution was in the (j) direction, until the field space 
path approached the valley bottom at (/) = 0. In cases were Wq is the dominant 
contribution to the energy density, this evolution was shown to lead to a positive 
spike in /nl with its maximum scaling as^"^ 

/NLlMax^O.Se/ -— , (14) 

A* 

which was verified for a concrete model. ^'^ Again we clearly see the fine tuning needed 
to generate a large /nl, and that for valleys a positive value of /nl is expected, a 
feature expected to apply to generic valleys. 

5. Models 

The discussions of the preceding sections become invaluable when we arc confronted 
with a concrete model of inflation. By identifying the ridge or the valley regions in 
the potential we can employ the above results to understand qualitatively how /nl 
evolves as these regions are traversed, and whether a 'natural' adiabatic limit will 
arise at the end of the evolution. Moreover, considering the conditions needed for 
a large non-Gaussianity at the adiabatic limit, and the estimates for a transitory 
large non-Gaussianity, we can identify initial conditions that give rise to appreciable 
non-Gaussianities. On the other hand, these analytic arguments can only inform us 
so far. In the literature there are examples of models for which analytic arguments 
suggest a large non-Gaussianity once a natural focusing region is reached, but where 
this only occurs after the slow-roll approximation ceases to accurately describe the 
dynamics. Moreover, there are examples of models for which no focusing region is 
available, but for which a large non-Gaussianity is possible as inflation ends. The 
various possibilities were classified in the introduction. To study these cases fully 
numerical simulations are essential. 

In our recent work^"^ we performed simulations of a variety of models, using a 
numerical implementation of the 6N formalism, confirming the usefulness of the 
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analytic estimates we have developed, but also highlighting how sensitive /nl can 
be to the exact time an adiabatic limit is reached, i.e. before or after slow-roll ends. 
In models for which no focusing region exists, we also showed a strong sensitivity 
of /nl on the time scale of reheating. 

Here we introduce an additional model, which highlights the need for numer- 
ical simulations, and also provides a new example of a model with a large non- 
Gaussianity after natural focusing. In this case the model can produce either a 
positive or negative asymptotic value of /nl depending on the initial conditions. 
We take a sum-separable potential of the form 

W{c^, X)^W„ + ]^ml4>^ + gx + ^Xx" + |Mx' (15) 

where W^, "t,^, t/, A and ^ are all positive constants. There is an inflection point 
in V{x) at X = 0. Wq and ^ are fixed by the requirement that there is a minimum 
with W^ = 0, at Xmin = — ''' where r is taken to be positive. These requirements 
impose two relations between the model parameters. This model clearly belongs to 
the class of potentials which contain a focusing region. 

An explicit example is given by g = (10~'*/0.18^)?7i|Mpi, A = 
(100/0. 18^)m?,/Mpi and r = 0.14Mpi, with the value of m^ implicitly chosen such 
that the value of C is normalised to be compatible with CMB constraints^^. Fig. 1 
illustrates the evolution of /nl with the initial condition (jf ~ 16Afpi, which sup- 
ports close to 60e-folds of inflation, and x* = 0.00157\/pi which represents an initial 
field position just above the infiection point. The ridge-like shape of the inflection 
point results in a negative spike in /nl- Subsequently the trajectories converge into 
a valley, while slow-roll is still maintained, resulting in /nl temporarily becoming 
positive before eventually tending to its limiting value. Recalling the discussion of 
§3, and bearing in mind the initial conditions, we find that N^^^ = V^/V' is the 
dominant first derivative of N ^ and moreover V!! /V^ is large and positive (since the 
field is initially above the infiection point). We therefore expect a large negative 
asymptotic value of /nl, which is borne out by the full numerical evolution. The 
numerical and HCA limiting values are both shown in Fig. 1, which show excellent 
agreement, since slow- roll is maintained throughout the evolution. We note that 
placing the field initially at x* = — O.OOlSMpi, where V^ is negative, leads to nearly 
identical evolution, except that a positive /nl is ultimately reached (with the nu- 
merically calculated value of /nl — 9.9 in good agreement with the analytic HCA 
value of /nl = 10.4). Finally, choosing initial conditions very close to x* — 0> leads 
to a negligible /nl, which again is expected since V^ ^ Q. 

Now let us change the parameters so that instead of reaching the adiabatic limit 
before infiation ends, it is reached just afterwards. We choose g = 10~"'TO^Mpi, 
A = 100m3,/A/pi and r = 0.14Afpi, and the same initial conditions (0* = 16A/pi, x* — 
O.OOlSAfpi). It is clear that the analytic HCA value for /nl can no longer be trusted, 
but one might hope that it gives at least an indication of the true asymptotic value, 
particularly as the x field will source a final fraction of an e-fold of inflation as it 
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Fig. 1. Evolution of /nl for the potential (15) with parameter values and initial conditions given 
in the text. The solid line shows the full numerical evolution while the dashed line gives the analytic 
solution. The thin vertical line at Af fa 64 represents the analytically calculated end of inflation 
when slow roll completely breaks down. The fact that the adiabatic limit value of /jml = —10.1 
(seen in the insert) has been reached long before this point explains why the 'HCA' analytic value 
of /nl — 10.8 is a good approximation. 

rolls. Studying Fig. 2, this is seen not to be the case, with the analytic estimate 
(unchanged from the above case) being extremely inaccurate. 



6. Conclusions 

We have studied the evolution of non-Gaussianity in multiple- field inflationary mod- 
els, using analytical and numerical methods. We have shown that the descent of 
fields from a ridge or their convergence into a valley can result in significant growth 
of non-Gaussianity, with the two cases being distinguished by the sign of /nl- 

To concretely predict non-Gaussianities one can employ analytical expressions 
or numerical methods. Currently, however, the former rely on the slow-roll approx- 
imation which limits their applicability. In such cases we have demonstrated that 
numerical methods can be invaluable. 

To calculate observationally relevant values of parameters such as /nl, in prac- 
tice it is necessary that an adiabatic limit is reached where these parameters become 
constant. In order to determine the possible behaviours of /nl as it reaches this 
limit, and the techniques needed to calculate it accurately there, we have found it 
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Fig. 2. Evolution of /nl for the potential (15) with parameters and initial conditions given in 
the text. As before, the solid line shows the full numerical evolution, the dashed line the slow-roll 
analytic evolution, and the thin vertical line shows the analytically calculated time when slow- 
roll breaks down. A discrepancy between the analytic and numerical evolution is not unexpected, 
though its magnitude is perhaps a surprise. 

is useful to classify inflationary models according to whether the adiabatic limit is 
reached naturally by the convergence of field space trajectories during the slow-roll 
regime, naturally after slow-roll ends, abruptly, or where no focusing region exists, 
only after reheating occurs. We have summarised the recent results from our paper^^ 
concerning these classes of models and given a new illustration here using a new 
model. 

Finally, we note that the new example included in this work indicates that it is 
easy to construct two-field sum-separable models which exhibit large /nl of positive 
or negative sign, even when the adiabatic limit is reached naturally during slow-roll 
inflation. We note, however, that all such models exhibiting a large non-Gaussianity 
at the adiabatic limit appear to require a large degree of fine tuning in their initial 
conditions. 
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